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ABSTRACT 



Type I planetary nebulae (PNe) are defined as those with high He and N 
abundances (Peimbert & Torres- Peimbert 1983). These objects present in gen- 
eral bipolar geometries and have high stellar temperatures (Corradi & Schwarz 
1995, Torres-Peimbert & Peimbert 1997). In this paper we analyse the empirical 
methods for abundance determination in order to check if the He and N over- 
abundances in Type I PNe are a consequence of a geometrical effect, due to the 
bipolarity, or the ionization stratification, due to the stellar temperature. For 
this, we obtain simulated spherically symmetrical as well as bipolar nebulae, us- 
ing a 3D photoionization code. From the projected emission line intensities for: 
a) the whole nebula; b) for a slit crossing the nebula; as well as c) for different 
positions in the nebula, we applied the formulae used in the literature to obtain 
empirical abundances. These empirical abundances are then compared with the 
adopted ones. We show that empirical abundances depend on the particular line 
of sight covered by the observation and can simulate an overabundance and/or 
the presence of abundance gradients of He and N in planetary nebulae with high 
stellar temperature. The geometrical effects are also discussed. Systematic errors 
in abundance determinations by empirical methods are higher for the N/H ratio 
than for N/0. Thus, it seems better to use the N/0 value when discussing N rich 
objects. 
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1. Introduction 

The knowledge of chemical abundances in Planetary Nebulae (PNe) is important since 
they can be used as constraints for stellar evolutionary models, as well as for studies of 
the chemical evolution of galaxies. In fact, PNe abundances reflect the products of stellar 
evolution, as well as the enrichment of the interstellar medium. Chemical abundances in PNe 
are usually obtained using empirical methods (as for example in Peimbert & Torres-Peimbert 
1987, Kingsburgh & Barlow 1994). Empirical methods are commonly used since they can 
easily be applied to a great number of objects, from a small number of bright emission lines. 
Photoionization models also provide the gas chemical composition when applied to specific 
nebulae (for example, Harrington et al. 1982; Clegg et al. 1987; Gruenwald, Viegas & 
Broguiere 1997). 

Following their abundance, PN are classified in types (Peimbert 1978; Peimbert & 
Torres-Peimbert 1983; Faundez-Abans & Maciel 1987). Type I PNe are those with high 
He and N abundances, e.g. He/H > 0.125 and log(N/0) > -0.30 (Peimbert & Torres- 
Peimbert 1983). In a study of morphological and physical properties of PNe, Corradi & 
Schwarz (1995) state that bipolar nebulae have the hottest stars among PNe, and, except 
for two objects, all bipolar PNe for which chemical data are available are Type I. Another 
common characteristic of bipolar nebulae is that the [Nil] line intensity is stronger than Ha 
all over the nebula, apart from the very central region (Corradi & Schwarz 1993b). This 
characteristic is usually explained by a N overabundance or by shock heating (Corradi & 
Schwarz 1993a; Corradi & Schwarz 1993c). 

Predictions of stellar evolutionary models and the high values obtained for He and N 
abundances in Type I/bipolar planetaries have suggested that their central stars are more 
massive than other PN nuclei (Peimbert & Serrano 1980; Calvet & Peimbert 1983; Torres- 
Peimbert 1984; Corradi & Schwarz 1995). However, new evolutionary models fail to explain 
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the high abundances of He and N in Type I PNe (Marigo, Bressan & Chiosi 1996), although 
models including the "hot botton burning" effect and accounting for the chemical evolution 
on the AGB may lead to higher N/0 abundances (van den Hoek & Groenewegen 1997). 

Prom the characteristics of Type I/bipolar planetary nebulae, it can be noted that 
objects with the highest He and N empirical abundances are those presenting a bipolar 
geometry and having the highest stellar temperatures. In fact, objects with higher stellar 
temperatures (T*) tend to present higher He and N empirical abundances (see Pig. 1). In 
order to homogenize the data in the figure, only objects with abundances obtained from line 
intensities integrated in a slit across the nebulae are plotted. As it will be shown in the 
following sections, empirical abundance determinations depend on the line of sight covered 
by the observation. Empirical abundances for the objects plotted in Pig. 1 are from de 
Preitas Pacheco et al. (1992), Kingsburgh & Barlow (1994), and Costa et al. (1996). The 
stellar temperature corresponds to the Zanstra temperature of Hell. Values for T* are taken 
from the literature (Kaler 1983; Preite-Martinez & Pottasch 1983; Pottasch 1984; Shaw & 
Kaler 1985; Gathier & Pottasch 1988; Gathier & Pottasch 1989; Glcizcs, Acker, & Stenholm 
1989; Kaler & Jacoby 1989; Shaw & Kaler 1989; Kaler, Shaw, & Kwitter 1990; Stanghellini, 
Corradi, & Schwarz 1993). 

In this paper we analyze the effect of the geometry and of the stellar temperature on the 
empirical abundance determinations, in order to verify if the high abundances of He and N 
obtained for Type I PNe are real. In fact, several studies show that the stellar temperature 
and/or the geometry can affect the abundances obtained from empirical methods. As dis- 
cussed by Peimbert and collaborators (Peimbert 1967; Peimbert & Costero 1969; Peimbert, 
Luridiana, & Torres- Peimbert 1995), the gas temperature obtained from observed forbidden 
line ratios and used in abundance determinations must be corrected by temperature fiuctu- 
ations inside the nebulae. In fact, for typical conditions of PNe, temperature fluctuations 
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can be important when the ionizing central star is hotter than 10^ K (Gruenwald & Viegas 
1995). Furthermore, a study on spherically symmetrical homogeneous HII regions show that 
temperature fluctuations obtained from projected data vary with the line of sight to the neb- 
ula; the variation depending on the characteristics of the star and the nebula (Gruenwald & 
Viegas 1992). With a three-dimensional (3D) self-consistent photoioinization code, a model 
for the bipolar PN IC 4406 has been obtained, showing that its abundance does not fulfill 
the criteria for Type I, as this nebula is usually classified (Gruenwald et al. 1997). Such an 
analysis shows that specific and more realistic models which take into account the geometry 
of the object can lead to abundance results that differ from those obtained with empirical 
methods. 

In order to check for the effects of geometry and stellar temperature on the empirical 
method to obtain PNe chemical abundances, we apply our self- consistent 3D photoionization 
code to planetary nebulae. For a wide range of nebular and stellar characteristics, we obtain 
the physical conditions in each point of the nebula, and calculate the resulting line intensity 
ratios at different lines of sight to the nebula. From the calculated line intensity ratios we 
determine the abundances applying the empirical methods used in the literature. These 
abundances are then compared to those assumed in the models. The models are described 
in §2. The results are presented in §3 and the conclusions are outhned in §4. 

2. Empirical Abundances from Simulated Nebulae 

The 3D code has been described in a previous paper (Gruenwald et al. 1997). As in 
classical uni-dimensional photoionization codes, the input parameters are the stellar temper- 
ature and luminosity, as well as the nebular elemental abundances and density distribution. 
Typical parameters for planetary nebulae are assumed: stellar temperatures (T*) ranging 
from 2x10^ to 3x10^ K, and stellar luminosities (L*) from 3x10^ L© to 2x10^ L©. The 
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models are obtained for two sets of chemical abundances (uniform throughout the nebula): 
solar abundances (Grevesse & Anders 1989), Zq, and an average abundance for planetary 
nebulae (Kingsburgh and Barlow 1994), Ztyp, which includes He, C, N, O, Ne, S, and Ar; 
for the remaining elements (Mg, Si, CI, and Fe), we assume a value equal to a hundredth 
solar in order to take into account their presence in grains (Stasinska & Tylenda 1986). For 
spherically symmetrical nebulae, a constant density, nj^, ranging from 10^ to 10® cm^^ is 
assumed. In order to analyze PNe with bipolar apparence, an equatorial torus of denser 
gas close to the star is included. The presence of a dense equatorial region can explain the 
observed geometry of bipolar PNe (Calvet & Peimbert 1983; Icke, Preston, & Bahck 1989; 
Sahai et al. 1991; Corradi & Schwarz 1995; Gruenwald et al. 1997). The torus has the 
following characteristics: i) external radius, r^xti equals to R/10, where R is the maximum 
dimension of the ionized region; ii) internal radius, rj„t, equal to rext/lO; iii) height, h^isfc, 
equal to lint'-, iv) ntorus= 10 hh- The maximum dimension R corresponds to the distance 
from the star where the fractional ionization of H, H+/H, equals to 5%, along the torus axis. 
The torus is centered at the central star. These values for the torus dimensions and density 
are chosen just to analyze the geometry effect on the abundance determinations, and are not 
intended to represent all the bipolar nebulae. 

The code provides the emissivity of the hues in each point of the nebula. Line intensity 
ratios can then be obtained: a) for different positions in the projected nebula; b) for the 
whole nebula; c) for a given slit across the nebula. 

We assume that the calculated line intensity ratios simulate observed ratios. Therefore, 
chemical abundances can be obtained by empirical methods using these line ratios. Com- 
paring the calculated abundances with those adopted for the simulated nebula, the possible 
effects of geometry and stellar temperature on the empirical abundance calculation can be 
established. The calculations are performed for different sets of emission line intensity ratios 
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corresponding to cases (a), (b), and (c), listed above. 

In order to obtain empirical elemental abundances, the ionic fractional abundance must 
be computed. For ions producing detected emission lines, the ion abundance relative to H+ is 
obtained from the observed line intensity relative to H/3 and the corresponding line emissivity 
ratio. Line emissivities depend on the gas density and temperature. The electronic density 
is obtained from the [SII] A6718/A6733 line intensity ratio, since this ratio depends weakly 
on the gas temperature (Osterbrock, 1989). Knowing the density, the gas temperature for 
the high ionization region {T[oiii]) can be obtained from the [0111] AA4959+5007]/A4363 
hne ratio, while for the low ionization region (T[jv//]) from the [Nil] AA6548+84]/A5755 
line ratio. Line emissivities, electronic density, and temperature are calculated assuming a 
five-level atom. 

A correction must be done for ions which are present in the gas and do not produce 
detectable emission lines. Such a correction can be done using the corresponding "ionization 
correction factor" (icf) for each element. We apply the zc/from Peimbert & Torres- Peimbert 
(1987), as generally used in the literature. Empirical abundances for He, O, and N are then 
obtained from the following equations: 



He 
If 



He+ + He++ 
H+ 



(1) 




(2) 




(3) 



where 



icfiO) 



He+ + He++ 



(4) 



-8- 



The ionization correction factor for O, icf{0), accounts for ionization stages higher than 
O"^"*". In equation (3), the last two terms correspond to the correction due, respectively, to 
N^"*" and to higher stages of N. 

The abundance of He"*" and He''"+, relative to H+, are determined respectively from the 
Hel A5876 and Hell A4686 line intensities relative to H/3. The ratio of line emissivities for 
these lines were calculated following Barker (1978), using the values obtained for the gas den- 
sity and temperature. A correction due to coUisional excitation of He"*" is taken into account 
following Peimbert & Torres-Peimbert (1987). The 0+, 0++, and abundances, relative 
to H+, are obtained respectively from [Oil] AA3727+29/H/3, [OIII] AA4959+5007/H/3, and 
[Nil] AA6548+84/H/3 hne ratios. 

3. Results 

In the following, the emission line ratios from simulated nebulae and the corresponding 
empirical chemical abundances are discussed for cases (a), (b), and (c) listed above. 

3.1. Spherically symmetrical nebulae 

The behavior of the [Nil] AA 6548-|-84/Hq! line intensity ratio with the projected distance 
to the center of the nebula is shown in Fig. 2a, while Fig. 2b gives the corresponding N/H 
abundance ratio. The simulations are performed with uniformly distributed density equal 
to 10^ cm Z — 'Zifyp, and stellar luminosity — 3x10"^ Lq. The projected distance to the 
center of the nebula is given in units of the radius of the ionized region. The curves are labeled 
by the stellar temperature, in units of 10^ K. The [Nil] /Ho; line ratio depends on the position 
of the line of sight to the nebula since different lines of sight cross different proportions of 
[Nil] and Ha emitting regions (Fig. 2a). The behavior of the [NII]/Hq; line ratio for different 
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lines of sight depends on the stellar temperature. Higher stellar temperatures lead to larger 
recombination regions, increasing the H"*" emitting volume relative to N"*" emitting volume. 
Therefore, there is a shift of the peak of [NII]/Hq; towards smaller rproj/R- with increasing 
T*, as shown in Fig. 2a. The maximum value for [Nil] /Ha depends on the value of the 
nitrogen abundance used in the models. However, as seen in the figure, the measured line 
ratio can substantially vary even for a nebula with a homogeneous abundance. 

Details on the [NH]/Hq! distribution on the projected nebula depend also on the stellar 
luminosity, on the nebular density, as well as on the chemical abundance. For increasing 
values of the nebular density, the region where [NH]/Hq: is higher than one are relatively 
narrower with the peak shifted to higher Tproj- The peak is also shifted to higher Tproj for 
increasing L^,. Lower values for the iron abundance increases the [Nil] /Ha ratio since Fe is an 
important cooler in the region where [Nil] lines are produced. This effect is more significant 
for lines of sight far from the central regions. 

Following Corradi & Schwarz (1993b), [Nil] stronger than Ha appears to be a common 
characteristic of bipolar nebulae. Our results show that even for a nebula with solar N/H 
abundance (which is lower than the average N/H for PNe), the measured [Nil] /Ha hne ratio 
can be higher than one. Indeed, a simulated nebulae with solar abundances and high T* 
(typical of Type I planetary nebulae) present [Nil] /Ha higher than one all over the nebula. 
For NGC 2440, which is a PN with a high value for T*, the [Nil] line is brighter than Ha, and 
the line ratio increases from 3 to 30, from the central to the outer zones (Icke et al. 1989). 
In fact, using characteristic parameters for this object (T^ = 2x10^ K, = 5x10^ Lq, and 
uh — 6x10^ cm~^), and for Z = Ztyp, we verify that the [Nil] /Ha ratio can reach values of 
the order of 26. Notice that these values are given just as an illustration, since these results 
correspond to a spherical symmetrical nebula, while NGC 2440 is clearly irregular (see, for 
example, Schwarz, Corradi & Melnick 1992). 
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The variation of the empirical N abundance in different lines of sight to a simulated 
nebula, is shown in Fig. 2b. The empirical values of N/H are plotted relative to the adopted 
abundance ratio. The empirical value obtained for the abundance depends on the detailed 
ionization structure of the nebula. Notice that for PNe whose central stellar temperature 
is high (T^, > 10^ K), the obtained values for different lines of sight can significantly vary, 
and can be higher than the adopted one in a extended region, mimicking an unreal over- 
abundance, as well as an abundance gradient. Therefore, studies of abundance gradients of 
nitrogen across a nebula applying empirical methods to the observed line intensities (Chu 
et al. 1991; Guerrero, Stanghelhni, & Manchado 1995; Guerrero, Manchado, & Serra-Ricart 
1996; Guerrero 1997; Corradi et al. 1997) can be misleading. For T^ < 10^ K, the resulting 
abundances are not very different from the adopted one, except in a narrow region near 
the borders. For such stellar temperatures, the volume of the H" zone is larger than the 
N° one. Lines of sight with increasing distance from the center include a higher fraction 
of low ionized regions. Therefore, empirical methods which do not account for neither N° 
nor H° overestimates the N abundance and must not be apphed in this case. For models 
with different nfj and/or L*, the discrepancies from the adopted abundance are somewhat 
different. However, the major effect comes from the stellar temperature. 

The variation of the [Nil] /Ho; with T^, is shown in Fig. 3a for the whole nebula and for 
a long and narrow slit across its central part. The two upper curves correspond to simulated 
nebulae with Ztyp, while the lower ones to Zq. The long-dashed and solid lines represent line 
ratios for the entire nebula, while the dot-dashed and short-dashed curves correspond to line 
intensities integrated in a slit. The nebular density and the stellar luminosity are the same 
as in Fig. 2. As expected, [Nil] /Ha values are lower for a slit than for the whole nebula. In 
fact, in the latter case there is a higher contribution of low ionized regions, increasing the 
[Nil] contribution relative to Ha, for a given T*. 
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Results for the empirical N/H relative to the adopted abundance are shown in Fig. 3b, 
for the case Z = Zq. For Z = Ztyp, the results are similar, except for T < 40000 K because 
of the effect of Fe on the gas cooling. The curves in Fig. 3b show that the empirical N/H 
abundances are overestimated for > 8x10^, while for lower the N abundance is lower 
than the true value. This behavior is related to the ionization correction factors used in the 
empirical calculations (see equations [1] to [4]). In fact, for T* < 8x10"* K the contribution of 
N"*""*" is underestimated since the correction for this non-observed ion is based on O"'""'" which 
has a smaller emitting volume than (the volume of the O"*" zone being greater than the 
N"*" one due to charge transfer reactions) . The main effect for high T^, is that the He"*""*" zone 
contains the N"*""*" zone; so, when applying the ionization correction factor for N"*""*", this ion is 
taken into account twice (both from the and He++ corrections) , increasing the obtained 
empirical N/H. Results for N/H from emission line intensities obtained with a slit crossing 
the nebular central region are greater than those where the whole nebula is observed. This 
comes from the fact that with such a slit the observed emitting volume of high ionized zones 
relative to low ionized zones is larger than for the whole nebula. 

The behavior of the 0/H empirical abundances with T* is similar to that of N/H, 
although less affected by the geometry. For example, when compared to N, the O abun- 
dance becomes greater than the true value for higher stellar temperatures. This leads to 
discrepancies of N/0 which increase with T^ (Fig. 3b). 

Since He/H is one of the criteria for classifying Type I PNe, this ratio is shown in Fig. 
4. For He/H, empirical methods do not take into account W nor He" abundances. However, 
for T* < 4x10^ K, there is an increasing He° zone, larger than the H" one; when He° is not 
taken into account, lower abundances for He are then obtained. For T* > 10^ K, the H° zone 
is increasingly larger than the He° one due to the high energy photons which can ionize He° 
very far into the nebula. If this effect is not properly taken into account, larger He/H are 
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obtained. 

For increasing values of uh, the discrepancies of empirical values of N/H increase, while 
those of He/H decrease (for increasing densities the volume of He" and H° tend to be equal). 
For example, for tlh = 10'^ cm~^, and same L^, as in Fig. 3b, the empirical N/H ratio is 
higher by a factor more than 3 for > 2 x 10^ K. Regarding the stellar luminosity, the effect 
is practically neghgible for N, but the discrepancies of He/H increase for lower L*. 

3.2. Bipoleir nebulae 

The effect of the torus and of the resulting bipolarity on the [Nil] /Ha hne ratio and 
on the corresponding empirical N/H are shown respectively in Fig. 5a and 5b. The curves 
correspond to T* = 10^ K, = 3x10^ L0, and Z = Ztyp. Solid lines correspond to a 
spherically symmetrical nebula with a uniform density {uh — 10^ cm~^) throughout the 
nebula. Dashed and dot-dashed curves show the result for a nebula with the same chemical 
abundance, T*, L*, and overall density, but including a torus around the central star. The 
torus density is ntoms — ^Ouh, and its dimensions are as described in §2. The dashed curve 
corresponds to a long sht along the lobes, and the dot-dashed one to a slit perpendicular 
to the torus axis, both passing through the central star. For a bipolar nebula the [Nil] /Ha 
ratio as well as the resulting N/H vary with the position of the line of sight to the nebula, 
due to its geometry. However, the effect due to the ionization distribution inside the nebula 
determined by the stellar temperature is much stronger, when the observations are made 
along the lobes, as it is often the case. For a given set of T*, L*, and overall density tlh, 
the main differences between the results for a bipolar and a homogeneous nebula are related 
to the torus characteristics, i.e, the contrast of n^ with the torus density, as well as to the 
torus dimensions. For a slit along the maximum extent of the lobes, the main differences 
appear for lines of sight crossing the torus. The opening angle of the torus will be important 
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for slits which make an angle with the torus axis. Any intermediate slit position will show 
a variation between the dashed and dotted-dashed curves of Fig. 5a and 5b. The results 
presented in these figures are an illustration of the ambiguity on the determination of any 
abundance gradients using empirical methods on bipolar objects. 

4. Conclusions 

The results of the preceding section show that empirical methods for abundance de- 
termination can simulate an overabundance and/or the presence of abundance gradients of 
He and N in planetary nebulae with high T*, even when applied to spherically symmetrical 
and homogeneous objects. This is also true for N/0. Recall that Type I PNe, defined as 
those nebulae which present an overabundance of He and N, have the hottest stars among 
PNe, and are generally bipolar. As shown in §3, both characteristics can induce an apparent 
overabundance when empirical methods are used. 

Following our results (Figs. 3a and 3b), systematic errors in the abundance determi- 
nations by empirical methods are smaller for N/0. Thus, it is better to use the N/0 ratio 
instead of N/H when discussing N rich planetaries. For central stellar temperatures in the 
range 10^ to 3x10^ K the N/0 ratio can be overestimated by a factor of 1.5. On the other 
hand, the average N/0 value for Type I PNe is about 2.8 higher than Type II PNe (Peimbert 
1990). So, we expect that, even with a better N/0 determination. Type I PNe may still be 
N rich compared to Type II PNe. 

Inference about abundances on planetary nebulae are important tests to stellar evo- 
lutionary models. Several arguments indicate that Type I PNe come from more massive 
progenitors (Torres- Peimbert & Peimbert 1997). In particular, high masses for the central 
stars of Type I PNe are generally invoked to explain their high He and N abundances since 
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the higher N/0 ratio, the higher the original mass of the central star. Some theoretical 
evolutionary models fail to explain these overabundances (Marigo et al. 1996). However, 
models accounting for the chemical evolution on AGB and including the "hot bottom burn- 
ing" effect are consistent with PNe abundances (van den Hoek & Groenewegen 1997). Our 
results show that He and N overabundances in Type I PNe may be lower than indicated by 
the empirical methods. The correction of the empirical abundances by those systematic er- 
rors change the constraints to the evolutionary models and may lead to different conclusions 
on the importance of the different processes. 

Furthermore, from the discussion in the preceding section, empirical methods applied 
to PNe with high stellar temperatures can mimic abundance gradients. Thus abundance 
gradients and their relation to multiple ejection phenomena must also be reviewed. 

Let us recall that observed line intensities are integrated on the line of sight, and different 
lines of sight cross different fractions of regions differing by their ionization state. Since 
lines originate in different regions, the observed image of a nebula depends on choice of 
the emission line. As shown in the preceding section, abundances obtained from observed 
line intensities depend on the line of sight. Therefore, empirical abundances obtained from 
line intensities averaged over different points in a nebula, or conclusions based on average 
abundances for a given nebula, may be incorrect. 

The effects shown in the present paper must be taken as a warning when discussing 
chemical enrichment and/or abundance gradients in a given nebula. Any conclusion about 
abundances, in particular in nebulae with a high star temperature, even with simple geome- 
tries, must rely on detailed photoionization models which properly account for the integration 
of emission line intensities on a given line of sight. Classical unidimensional photoionization 
models are suitable for spherically symmetrical objects, whereas a 3D photoionization code 
(Gruenwald at. al 1997) must be used for other geometries. 
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Fig. 1. — Empirical abundance ratios from the literature (see text) versus the stellar tem- 
perature for N/H (a) and He/H (b). 

Fig. 2. — [Nil] /Ha line ratio (a) and the empirical N/H (b) relative to the adopted value 
as a function of the line of sight. The projected distance to the center of the nebula is 
given in units of the nebula ionized radius. Curves correspond to simulated spherically 
symmetrical nebulae with = 3000 Lq, uh — 10^ cm~^, and Ztyp, and are labeled by the 
stellar temperature in units of 10^ K. 

Fig. 3. — (a) [Nil] /Ha versus T^. The two upper curves correspond to Z^^p, while the lower 
ones to Zq. Long-dashed and solid lines represent line ratios for the whole nebula, while 
dot-dashed and short-dashed curves correspond to line intensities integrated in a slit. The 
values of and uh are the same as in Fig. 2. (b) empirical N/H and N/0 relative to the 
adopted value (Z©). 

Fig. 4. — same as Fig. 3b, for He/H 

Fig. 5. — [NII]/Hq; (a) and N/H (b) relative to the adopted ratio as a function of the 
line og sight, for T^, = 10^ K. Solid lines correspond to a spherically symmetrical nebula, 
while dashed and dot-dashed correspond to a bipolar nebula of same stellar temperature, 
respectively for a slit along the lobes, and a slit perpendicular to the torus axis. 
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